NASA Technical Memorandum 104648 


Spherical Harmonics Analysis of the 
ECMWF Global Wind Fields at the 
10-Meter Height Level During 1985: 

A Collection of Figures Illustrating Results 


Braulio V. Sanchez 
Goddard Space Flight Center 
Greenbelt, Maryland 

Masahiro Nishihama 
Hughes STX 
Greenbelt, Maryland 



National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 

1997 




NASA Technical Memorandum 104648 


Abstract 

Half-daily global wind speeds in the east- west (u) and north-south (v) directions at the 10-meter 
height level were obtained from the European Centre for Medium Range Weather Forecasts 
(ECMWF) data set of global analyses. The data set covered the period 1985 January to 1995 January. 

A spherical harmonic expansion to degree and order 50 was used to perform harmonic analysis of 
the east-west (u) and north-south (v) velocity field components. 

The resulting wind field is displayed, as well as the residual of the fit, at a particular time. The 
contribution of particular coefficients is shown. The time variability of the coefficients up to degree 
and order 3 is presented. Corresponding power spectrum plots are given. Time series analyses were 
applied also to the power associated with degrees 0-10; the results are included. 



Braulio V. Sanchez and Masahiro Nishihama 


Acknowledgments 

We express our appreciation to Dr. Nikolaos K. Pavlis of Hughes STX for providing the Fast Fourier 
Transform software. We give thanks to Dr. Andrew Y. Au of Hughes STX for facilitating access to 
the ECMWF data set. This work was supported by the Space Geodesy Branch at the Goddard Space 
Flight Center and the TOPEX/POSEIDON project. 


iv 



NASA Technical Memorandum 104648 


Table of Contents 

Introduction 1 

Global Wind Fields 3 

Figures 1-10 5 

Time Variability of Selected Spherical Harmonics Coefficients 17 

Figures 11-45 19 

References 57 

Appendix 1 59 

Appendix II 61 




NASA Technical Memorandum 104648 


Introduction 


The general circulation of the ocean is driven by various types of forces, including surface wind 
stress. The results presented in this publication were obtained as part of an investigation of the wind- 
driven variability of the ocean circulation. It is possible that these results might be of interest to other 
investigators, not only in the limited context of ocean circulation, but also from the larger perspective 
of geodynamics. Earth rotation, and space geodesy. With such possibility in mind, it was decided to 
document the results. However, the intention is only to present the results obtained in graphical form; 
no interpretation is attempted here. 

Half-daily global wind speeds in the east-west (u) and north-south (v) directions at the 10-meter 
height level were obtained from the European Centre for Medium Range Weather Forecasts 
(ECMWF) data set of global analyses. The data set covered the period 1985 January to 1995 January. 
The original global grid is 2.5° latitude by 2.5° longitude, and the data set used in this investigation 
has been reprocessed and interpolated onto a 2° latitude by 2.5° longitude grid, as described by 
Schubert et al. (1990). 

A spherical harmonic expansion to degree and order 50 was used to perform harmonic analysis of 
the east-west (u) and north-south (v) velocity field components. The formal validity of this procedure 
was established by Simmonds (1974). Formal developments are given in Appendix I. 

The software used in the analysis implements a technique developed by O. L. Colombo (1981), 
based on the Fast Fourier Transform. 
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Global Wind Fields 

Figures 1-10 show the vectors for the wind fields corresponding to various cases. With the exception 
of figures 7 and 8, the time epoch is January 1, 0 hour GMT, 1985. This is the time of the first data 
entry for 1985. 

Figure 1 shows the wind field as obtained from the ECMWF data files, the maximum magnitude is 
22 meters/second. This wind field was fitted by means of an expansion in spherical harmonics 
including coefficients up to degree and order 10. The resulting field is shown in figure 2; maximum 
wind speed is 12 meters/second. The difference between the data and the fit is displayed in figure 3; 
the largest discrepancy has a magnitude of 16 meters/second. 

Figure 4 exhibits the wind field obtained by fitting the data with spherical harmonics up to degree 
and order 50, note the improvement as compared to the degree- 10 fit; the maximum wind speed is 
now 20 meters/second. The residuals are shown in figure 5; the maximum error has a magnitude of 
4.5 meters/second. Note that the largest discrepancies are concentrated around the poles. 

Figures 6 and 7 exhibit the wind fields from the ECMWF data files at January 1, 0 hour GMT, 1985 
and at July 1, 0 hour GMT, 1985; the wind vectors are plotted with equal length in order to show 
clearly the wind direction, the magnitudes are indicated by the shading, the lower (higher) wind 
speeds corresponding to the darker (brighter) shades. Figure 8 depicts the average wind field for 
1985. 

Figures 9 and 10 are included to display the effect of some of the coefficients in the spherical har- 
monic expansion. Figure 9 brings to view the wind field corresponding to the C 10 coefficient by 
itself. Figure 10 makes apparent the wind field produced by the sum of the C„ and S n coefficients. 
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FIGURES 1-10 
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Figure 1. Wind Field at 1985, January 1, 00. As Obtained from the ECMWF Data Set, 
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Figure 2. Wind Field at 1985, January 1, 00. As Obtained from the Spherical Harmonic Expansion to Degree 
and Order 10. 
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Expansion to Degree and Order 10. 
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Figure 4. Wind Field at 1985, January 1, 00. As Obtained from the Spherical Harmonic Expansion to Degree 
and Order 50. 
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Figure 5. Vector Difference Between the ECMWF Wind Field, and the Wind Field from the Spherical Harmonic 
Expansion to Degree and Order 50. 
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Figure 6. Wind Field Directions at 1985, January 1, 00. As Obtained from the ECMWF Data Set. 
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Figure 7. Wind Field Directions at 1985, July 1, 00. As Obtained from the ECMWF Data Set. 
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Figure 8. Average Wind Field Directions for 1985. As Obtained from the ECMWF Data Set. 
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Figure 9. Wind Field at 1985, January 1, 00. As Obtained from the Spherical Harmonic Expansion, Using Only 
the C, n Coefficient. 
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Figure 10. Wind Field at 1985, January 1, 00. As Obtained from the Spherical Harmonic Expansion, Using 
Only the C u and S u Coefficients. 
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Time Variability of Selected Spherical Harmonic Coefficients 

Figures 1 1-26 portray the time variability of the spherical harmonic coefficients up to degree and 
order 3. Each figure shows the results for the east-west (u) and north-south (v) velocity components. 
There are 730 values over the course of a year for each spherical harmonic coefficient. These time 
series were analyzed by means of Fourier expansions in sine and cosine terms, as presented in 
Appendix II. Plots of power as a function of frequency are given in figures 27-34. In general, most 
of the power is concentrated in the low frequency terms, especially the zero frequency component. 
The east-west velocity component predominates over the north-south, with the exception of the 
coefficients C^, C 30 , and S 3r 

The time series analysis was applied also to the power associated with degrees 0-10, as given by the 
last equation in Appendix I. The results are dispayed in figures 35-45. Each figure presents the 
results for the east-west (u) and north-south (v) velocity components. As shown in part (a) of the 
figures, most of the power is contained in the zero-frequency term. Part (b) makes evident the power 
in the other frequencies by setting the constant part equal to zero. 
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FIGURES 1 1-45 
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Figure 12b. The Coefficient C 10 , Year 1985. The North-South Velocity Component (v). 
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Figure 18a, The Coefficient C 22 , Year 1985. The East- West Velocity Component (u). 
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Figure 18b. The Coefficient C 22 , Year 1985. The North-South Velocity Component (v). 
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U10M COEFFICIENT S22 IN 1985 
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Figure 19a. The Coefficient S 22 , Year 1985. The East-West Velocity Component (u) 
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Figure 19b. The Coefficient S 22 , Year 1985. The North-South Velocity Component (v) 
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Figure 22b. The Coefficient S 31 , Year 1985. The North-South Velocity Component (v). 


























Figure 25b. The Coefficient C 33 , Year 1985. The North-South Velocity Component (v) 



















NASA Technical Memorandum 104648 

U/V 10M POWER SPECTRUM OF COO OVER 1985 



Figure 27a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C,,,,. 
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Figure 27b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C 10 . 
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U/V 10M POWER SPECTRUM OF C20 OVER 1985 
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Figure 28a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C^. 



Figure 28b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C^. 
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Figure 29a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C n . 



Figure 29b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S n . 
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Figure 30a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C 21 . 
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Figure 30b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S 21 . 
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Figure 31a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C 22 . 



Figure 31b. Power Spectrum. Frequencies: 0—50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S 22 . 
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Figure 32a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East-West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C M . 
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Figure 32b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S 31 . 
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Figure 33a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C 32 . 

U/V 10M POWER SPECTRUM OF S32 OVER 1985 
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Figure 33b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S 32 . 
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U/V 10M POWER SPECTRUM OF C33 OVER 1985 
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Figure 34a. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient C 33 . 



Figure 34b. Power Spectrum. Frequencies: 0-50. Year 1985. Solid Line: East- West Velocity 
Component (u). Dashed Line: North-South Velocity Component (v). The Coefficient S 33 . 
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DEGREE 0 POWER SPECTRUM FOR 1985 



Figure 35a. Degree 0, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 


DEGREE 0 POWER SPECTRUM FOR 1985 
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Figure 35b. Degree 0, Power Spectrum. Year 1985. Solid Line: East-West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 1 POWER SPECTRUM FOR 1985 
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Figure 36a. Degree 1, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 



Figure 36b. Degree 1, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 2 POWER SPECTRUM FOR 1985 
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Figure 37a. Degree 2, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 
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Figure 37b. Degree 2, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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Figure 38a. Degree 3, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 



Figure 38b. Degree 3, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 4 POWER SPECTRUM FOR 1985 
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Figure 39a. Degree 4, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 
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Figure 39b. Degree 4, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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Figure 40a. Degree 5, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 



Figure 40b. Degree 5, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 6 POWER SPECTRUM FOR 1985 
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Figure 41a. Degree 6, Power Spectrum. Year 1985. Solid Line: East-West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 


Figure 41b. Degree 6, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 7 POWER SPECTRUM FOR 1985 

'll. 

10 . 

9.8 

9.3 ■ 

8.7 
8.2 

7.6 

7.1 

6.5 

oc 6.0 ■! 

Lii 

§ 5.5 ■! 

O 

4.9 ■: 

4.4 ■■! 

3.8 • ; 

3.3 ; 

2.7 

2.2 

1.6 
1.1 
.55 

.0 Lb- — * * * 

0 10 20 30 40 50 

FREQUENCY 


Figure 42a. Degree 7, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 



Figure 42b. Degree 7, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 8 POWER SPECTRUM FOR 1985 
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Figure 43a. Degree 8, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 
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Figure 43b. Degree 8, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u) 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 













Figure 44a. Degree 9, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 



Figure 44b. Degree 9, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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DEGREE 10 POWER SPECTRUM FOR 1985 
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Figure 45a. Degree 10, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-50. 
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Figure 45b. Degree 10, Power Spectrum. Year 1985. Solid Line: East- West Velocity Component (u). 
Dashed Line: North-South Velocity Component (v). Frequencies: 0-40, Power Set to Zero at Zero 
Frequency. 
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APPENDIX I 

HARMONIC ANALYSIS OF FUNCTIONS DEFINED ON THE GLOBE 

Let f(9, X) represent a function defined globally at a set of discrete points, that is, 

f(e, x)=m i ,x.) 

i = 0, 1, ... , n-1 
j =0, 1, ..., m-1 

A. : longitude A. - A. = 2jt/m 

9. : colatitude 9. +) -9. = n/n 

Assume that f(9. A) can be expanded in a series of spherical harmonics, 
f(9. A) = I n I m \ (C nm cos(mA) + S nm sin(mA) } P nm (9) 

where, 

P (0) : fully normalized Legendre functions of degree n, 
order m, i.e. 

j J P nm (9) ^ cos(mA), sin(mA) ^ P k] (9) \ cos(lA), sin(lA) sin(9) d9 dA 

= W S nL 5 m. 

Then, 

C nm = (1/4 tc) J j f (9, A) P nm (9) cos(mA) sin9 d9 dA 

S nm = ( 1/471) j 1 f(9, A) P nm (9) sin(mA) sin9 d9 dA 

The power associated with the terms of degree n is given by, 

w =x (C ) 2 + (S ) 2 

n m v nm y v nm 7 
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APPENDIX II 

POWER SPECTRUM OF FUNCTIONS DEFINED ON THE GLOBE 

Let F denote a function such as C , S , or W , as defined in Appendix I. If certain basic 

nm nm n r r 

conditions are satisfied, then F can be expanded in a Fourier series: 

F = a 0 + a, cos(t) + b, sin (t) + ^ cos (2t) + b 2 sin (2t) + . . . + a 365 cos(365t) + b 365 sin(365t). 
for 0 < t < 2rc 

If the data covers a time span of one year, then the zero frequency term corresponds to the average 
over the year, frequency 1 corresponds to a periodicity of one year, frequency 2 corresponds to a 
periodicity of 6 months, and so on. The power corresponding to each frequency is given by: 

P i = (a j ) 2 + (b.) 2 
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